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ABSTRACT 

Deep X-ray imaging spectroscopy of the bright pulsar wind nebula 3C58 confirms the existence of 
an embedded thermal X-ray shell surrounding the pulsar PSR J0205-I-6449. Radially resolved spectra 
obtained with the XMM-Newton telescope are well-characterized by a power-law model with the 
addition of a soft thermal emission component in varying proportions. These fits reproduce the well- 
studied increase in the spectral index with radius attributed to synchrotron burn-off of high energy 
electrons. Most interestingly, a radially resolved thermal component is shown to map out a shell-like 
structure « 6' in diameter. The presence of a strong emission line corresponding to the Ne IX He-like 
transition requires an overabundance of ~ 3 x [Nc/NCq] in the Raymond- Smith plasma model. The 
best-fit temperature kT ^ 0.23 keV is essentially independent of radius for the derived column density 
of TVh = (4.2 ± 0.1) X 10^^ cm~^. Our result suggests that thermal shells can be obscured in the early 
evolution of a supernova remnant by non-thermal pulsar wind nebulae emission; the luminosity of 
the 3C58 shell is more than an order of magnitude below the upper limit on a similar shell in the 
Crab Nebula. We find the shell centroid to be offset from the pulsar location. If this neutron star 
has a velocity similar to that of the Crab pulsar, we derive an age of 3700 yr and a velocity vector 
aligned with the long axis of the PWN. The shell parameters and pulsar offset add to the accumulating 
evidence that 3C58 is not the remnant of the supernova of CE 1181. 

Subject headings: stars: individual (PSR J0205-I-6449, 3C58) — ISM: supernova remnant — stars: 
neutron — X-rays: stars — pulsars: general 



1. INTRODUCTION 

A long standing puzzle in supernova physics is the ap- 
parent lack of an associated thermal shell surrounding a 
few young pulsars with bright relativistic wind-powered 
nebulae (PWNe). Though relatively rare, with fewer 
than ten examples known, the lack of a supernova rem- 
nant shell is at odds with our current understanding of 
supernova evolution and neutron star formation. The 
canonical examples, the Crab and 3C58, are both pu- 
tative historical remnants less than 1000 yrs old. To 
date, no evidence has been found for thermal emission 
associated with the Crab Nebula pulsar despite repeated 
searches (Seward et al. 2006 and references therein). 
However, recent XMM-Newton and Chandra observa- 
tions have identified thermal emission within the 3C58 
remnant (Bocchino et al. 2001; Slane et al. 2004). The 
availability of a large quantity of XMM archival data 
on this source derived from calibration observations al- 
lows us to construct very sensitive images and spectra to 
explore the evidence for a surrounding SNR shell with 
unprecedented sensitivity. 

The morphology of 3C58 is well-documented in both 
the X-ray and radio bands (Slane et al. 2004 and Bi- 
etenholz et al. 2006 and references therein). It is center- 
filled with a 66 msec pulsar near the center, axisymmetric 
lobes, and an elaborate network of wisps and filaments. 
Failure to find a radio shell around 3C58 led Reynolds 
& AUer (1985) to conclude that there is no evidence for 
its interaction with an external medium. Coupled with 
the low velocities of the optical filaments (Fesen 1983) 
compared to the mean expansion velocity required for 
the remnant to reach its current extent in 820 yr, these 
results led to suggestions that the remnant was consider- 
ably older. Over the past twenty years, further evidence 



has accumulated that the remnant's age is inconsistent 
with an origin in SN1181. The images and spatially re- 
solved spectra we present below add incrementially to 
the case and provides a testable prediction concerning 
the pulsar velocity which could support or refute the his- 
torical association. 

In §2 we describe the available observations and discuss 
their analysis. Section 3 presents the detailed spatially 
resolved spectroscopy which allows us to detect clear ev- 
idence for a shell of thermal emission, overabundant in 
Neon, and extending to the edge of the synchrotron neb- 
ula; we interpret this result in §4, concluding that the 
pulsar is offset from the center of the symmetrically ex- 
panding shell implying a velocity aligned with the sym- 
metry axis of the remnant. If the pulsar has a Crab-like 
velocity, this offset provides another argument against an 
association with SN1181. 

2. OBSERVATIONS AND RESULTS 
2.1. The observations 

A total of ten pointed observations of 3C58 are avail- 
able in the public archive of data obtained with the 
XMM-Newton X-ray Observatory (Jansen et al. 2001). 
Of these, a set of eight pointings, acquired on 2002 
September 11 — 13 to perform an in-orbit calibration 
of the telescope mirror vignetting function (Lumb et al. 
2004), used imaging modes with a sufficiently large field- 
of-view to cover fully the 3C58 nebula. The duration of 
these observations varied between 17 — 33 ks. An obser- 
vation log is presented in Table 1. 

In this study we concentrate on data obtained with the 
European Photon Imaging Camera (EPIC; Turner et al. 
2001) which consists of three CCD cameras, the EPIC-pn 
and the two EPIC-MOS imagers. The EPIC-MOS cam- 
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TABLE 1 

XMM-Newton Observation log for 3C58" 



ObsID 


Date 


Pointing 


pointing 


Exposure 


Exposure 






R.A 


Decl. 


EPIC-MOS* 


EPIC-pn 




(UT) 


(J2000) 


(J2000) 


(ks) 


(ks) 



153752201 


2002-09-11 04:29:35 


2h4m43 gs 


-1-64° 51' 13^2 


36777 


15579 


153751801 


2002-09-11 13:09:58 


2''5™03.0^ 


-1-64° 47' 38^6 


39561 


16328 


153752501 


2002-09-11 20:03:41 


2''5™23.4= 


-|-64°43'52'.'3 


35195 


14721 


153752401 


2002-09-12 03:20:44 


2''6™32.4= 


-t-64°48'06'.'8 


26184 




153752101 


2002-09-13 04:22:11 


2'\5™38.0= 


-t-64°49'40'.'0 


17363 


3022 


153751701 


2002-09-13 10:57:39 


2ii5™18.9= 


-|-64°53'23y7 


13515 


6042 


153751901 


2002-09-13 17:51:22 


2^5^13.1= 


+64°51'41^'8 


34114 


14569 


153752001 


2002-09-13 23:45:05 


2''5™52.6= 


-l-64°55'2r'7 


48631 


21038 



^AII observations were obtained in pointed imaging mode with the thin transmission window in the focal plane; the EMOS and EPN data 
were collected with the PrimeFullWindow and PrimeLargeWindow submodes, respectively 
'^EPIC exposure times are after filtering; the EPIC-MOS exposure times are for the sum of the two MOS cameras. 




R. A. (JZOOO) R. A. (JZOOO) R. A. (JZOOO) 

Fig. 1. — Narrow-band X-ray images of 3C58 acquired with the JfMM-Newton EPIC cameras. The energy range for these exposure- 
corrected images separates out three distinct sources of emission, Left) the unresolved point- like pulsar emission surrounded by the highest- 
energy and most intense portion of the PWN (7.0 — 8.0 keV), Middle) the pulsar wind nebula in the medium band (2.0 — 3.0 keV), and 
Right) the symmetric shell-like structure in the soft X-ray band (0.5 — 1.0 keV). The pulsar wind nebula image also includes pulsar emission 
whereas the shell-like structure includes both pulsar and PWN emission components. In each image, the intensity is cropped at 0.6 of the 
peaJj intensity and plotted with square-root scaling to highlight the fainter emission. 



eras were operated in FullFrameMode mode, in which 
the full 30' diameter field-of-view is read out every 2.7 s. 
The EPIC-pn was operated in PrimeLargeWindow mode, 
which provides 78 ms time resolution over a 6f9 x 13f5 
field-of-view. For all observation, the thin transmission 
filter was placed in the focal plane. The EPIC instru- 
ments are sensitive to X-rays in the 0.2 — 10 keV energy 
range. The point spread function (PSF) of the mirror 
modules has a FWHM of 6 — 15" depending on energy 
and is fully oversampled by the CCD pixel size in each 
instrument. 

2.2. Basic analysis 

We analyzed data produced by the standard processing 
(SAS version 20020507.1701). The photon event lists 
were initially filtered for periods of high background, typ- 
ically fiare events corresponding to enhanced solar wind 
activity. To identify these flares, an iterative clipping 
method was used to generate the acceptable time inter- 
vals using the following method. We first produced a 
lightcurve histogram in 100 s steps for the whole instru- 



ment and calculated the mean. Intervals with an en- 
hanced count rate > 3cr above the mean were discarded 

and a new lightcurve was generated. This process was 
repeated until all points above the converging threshold 
were removed. The good time intervals (GTIa) thus gen- 
erated were used to filter the revised event file. A total 
of ~ 343 ks of good exposure time was acquired from the 
collective data set; the net exposure times for each instru- 
ment are listed in Table 1. Finally, in creating spectra 
and images, we used the standard SAS screening flags 
and selected only CCD photons with PATTERN < 4 and 
< 12 for the EPIC-pn and EPIC-MOS, respectively Bad 
CCD pixels and columns were taken into account. We 
verified that spectra were not affected by photon pile-up 
in any of the observations. 

2.3. Spatial analysis 

To search for evidence of a faint supernova rem- 
nant surrounding the pulsar, we generated exposure- 
corrected, narrow-band images including data from all 
three instruments in order to obtain the deepest possible 
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Fig. 2. — Evidence of a thermal shell surrounding the pulsar 
in 3C58. Normalized background-subtracted radial profiles for five 
adjacent energy bands from 0.3 — 8.0 keV obtained with the XMM- 
Newton EPIC CCD detectors. A prominent plateau of emission is 
seen for the lowest energy profile. A shell-like structure {dashed- 
line) is found for the radial profile of the 0.3 — 1.2 keV band image 
after subtracting off the normalized emission in the next highest 
band. 



0.05 




Powei — law + 
Thermal Model 



t^*t 

1 1 1 — h- 1— H 



^+-*^+■H;^*M^tfpJ%lt^^ 



Powet — law + J 
Thermal + Ne IX Model 



ENERGY (KEV) 

FlK. 3. Evidence for a tlicriual component with an overabun- 
dance of Neon in 3C58. The solid line in the upper panel represents 
a best-fit power-law plus Raymond-Smith plasma emission model 
to the XMM-Newton EPIC-MOS spectra (crosses) from annular 
regions (2.'0 — 2f5) centered on the pulsar in 3C58. The dotted 
line shows the power-law component only. Residuals from three 
different fits are displayed in the lower three panels. The first of 
these shows a strong excess above the pure power-law model below 
1.2 keV. The next panel down shows the residuals from a model 
that includes an equilibrium thermal R-S plasma with solar abun- 
dances. The strong line remaining at ~ 0.9 keV corresponds to 
a complex of emission lines from Helium-like Neon (Ne IX). The 
bottom panel shows a fit in which the Neon abundance is three 
times solar. Model parameters are given in Table 2. 



image. Figure 1 compares the image of 3C58 in three 
energy bands to highlight the three spatial components 
which are found to make up the composite morphology: 
a shell-like nebula (0.5 — 1.0 keV), the pulsar wind neb- 
ula (2.0 — 3.0 keV), and the point-like pulsar emission 



(7.0 — 8.0 kcV). This is best illustrated by comparing the 
radial profiles in the five narrow energy-bands shown in 
Figure 2. In addition to the steep falloff of flux with ra- 
dius, it is again clear that the extent of the emission from 
3C58 is strongly dependent on energy-band. In particu- 
lar, the radial profile in the softest energy band displays 
a largo departure from the higher energy profiles between 
2' and 5' from the remnant center. Subtracting emission 
in the 1.2 — 2.0 keV band from the lowest energy image 
shows that the soft emission has a radial profile consis- 
tent with that of a limb-brightened shell peaking ~ 3' 
from the remnant center. 

To investigate further the nature of this shell-like emis- 
sion we generated spectra from concentric annuli 30" 
wide to search for radial variations. These annuli were 
centered on the location of the pulsar, the presumed ori- 
gin of the wind and likely center of any SNR expan- 
sion. To estimate the background in each annulus, we 
collected photons from within a region with ~ 3f5 < 
r < 4f0. Although the background is somewhat distant 
from the inner source regions, the high surface brightness 
of the source makes background subtraction less impor- 
tant there. We also extracted photons from the eastern 
"lobe" region, using a circle of radius 75" centered at 
coordinates RA = 02:06:14.65 and Dec = 64:49:31.9; the 
background for this region was determined from a pair of 
apertures straddling the source region to the north and 
south. 

2.4. Spectral analysis 

Spectra from each observation were extracted from 
each annular region using standard channel binning for 
each instrument, and a set of response matrices were gen- 
erate for each ring using the standard prescription for dif- 
fuse emission. The spectra and response matrices were 
then summed over all observations for each instrument 
following the method used in addascaspec^. The two 
sets of MOS spectra were then combined. Finally, all 
spectra were grouped into bins containing a minimum 
of 400 counts and fitted using the xspec spectral fitting 
package. 

To fit the summed spectra in each annulus we used 
the following procedure. In all cases, the EPIC-pn and 
EPIC-MOS spectra were fitted simultaneously, with the 
model normalizations allowed to be independent to ac- 
count for remaining flux calibration differences ( ^ 10%) 
between the two instruments (Saxton 2003). The spectra 
were initially fitted in the 1.5 — 8.0 keV energy band with 
a power-law model for the non-thermal emission from the 
pulsar wind nebula. This model included the effects of 
interstellar absorption, although this is not an important 
component in the fitted band. The best-fit values for the 
power-law index and flux are recorded in Table 2. As 
illustrated in Figure 3, the spectrum above 1.2 keV is 
well-characterized by a power-law model. However, at 
the lower energies, we require additional emission com- 
ponents for the larger annuli. We account for this excess 
by adding in a Raymond-Smith thermal plasma emission 
component to the model and extending the fitting range 
to 0.5 — 8.0 keV. In fitting this new model, the power- 
law indices were fixed to the values derived at high ener- 

^ addascaspec is part of the FTOOLS software package available 
at heasarc . gsf c. nasa . gov/docs/sof tware . html 
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Fig. 4. — Left Panel: Residual from the best fit to the XMM-Newton EPIC-MOS spectrum of 3C58 in units of sigma as a function of 
radius. Data are accumulated in concentric aiinuli centered on the pulsar as defined in Table 2 and are ordered by decreasing radius (top 
to bottom). Each plot shows data from the combined EPIC-MOS spectrum from eight observations of the pulsar; a similar result is found 
for the EPIC-PN camera (not shown). The data in all cases are well-fit by a two-component power-law plus Raymond-Smith model with 
the parameters given in Table 2. Right Panel: As for adjacent panel, but with the thermal component normalization set to zero to show 
this contribution as a function of radius. The flux of the thermal component as a function of radius is also displayed in Figure 6. 
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Fig. 5. XMM-Ncwton EPIC spectrum of the 3C58 lobe region 
(see text) fitted to a power-law model. The bottom panel displays 
the residuals from the best fit model with P = 2.88 ± 0.05 keV. 



gies. As shown in the third panel of Figure 3, however, 
a strong Une feature at ~ 0.9 keV remained unmodeled. 
This emission hue corresponds to the He-hke Ne IX tran- 
sition and requires an overabundance ~ 3 x [Ne/Ne^] in 
the Raymond-Smith plasma model throughout the neb- 
ula. 

Our final model includes a Gaussian line component to 
allow for the Neon emission contribution. For these fits 
the Nff was fixed to the value derived from the central 
annulus, the region with the most counts and in which 
the non-thermal emission strongly dominates over any 
thermal emission. In all cases, the final spectral fits pro- 
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Fig. 6. — Spectroscopic detection of a thermal shell surround- 
ing 3C58. Shown are the normalized radial profiles of the com- 
ponent flux from a power-law plus variable abundance Raymond- 
Smith plasma emission model. The observed fluxes are given in the 
0.5-10.0 keV energy band using the EPIC-MOS cameras on-board 
XMM-Newton; a similar result is found using the EPIC-pn cam- 
eras (see Table 2). Notice the change of scale by a factor of ~ 500 
between the ordinate axes; less than 4% of the observed remnant 
flux is in the thermal component, explaining why it has proven so 
illusive. 



duce an acceptable fit statistic, with the exception of the 
region between r = 2 — 2.5', and perhaps the next interior 
one. This can be attributed to unmodeled structure be- 
low 0.8 keV, likely additional weak thermal line emission 
which we lack the sensitivity to model with the current 
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Fig. 7. — Photon index for the best fit power-law model, 
Raymond-Smith plasma temperature, and Ne abundance as a func- 
tion of radius. The power-law index increases monotonically from 
r = 1.9 to r = 2.8, consistent with previous reports in Torii et al. 
(2000) and Slanc ct al. (2004). The index derived from emission 
in the eastern lobe up to 4' from the center is shown on the right 
axis. Both the plasma temerature and the Neon abundance show 
no significant radial dependence between l.'O and 3.'2 



data. Tabic 2 gives the best fit values for the tempera- 
ture and flux in each annulus after fixing the F and N//, 
while Figures 4 demonstrates the dramatic change in the 
thermal contribution as a function of radius and the high 
quality of the fits in the composite model. 

The overall profile of the remnant is asymmetrical, 
with a lobe extending to the east. In order to deter- 
mine whether this is emission is powered strictly by an 
asymmetric relativistic wind from the pulsar (cf. Vela, 
Helfand et al 2001; B1509-58, Gaensler et al. 2002) 
or represents the breakout of the hot gas into a re- 
gion of low surrounding density as is seen in some shell- 
type remnants, we analyzed the spectrum of the emis- 
sion extracted from the lobe region defined above. Fig- 
ure 5 displays the spectrum along with residuals from a 
power law fit which fixed the column density at the value 
used above for the remnant as a whole. The spectrum 
is well-characterized by a steep power law with index 
r = 2.88 ± 0.06 with no hint of any thermal emission at 
low energies. The spectral index is the steepest value we 
measured, consistent with its location farthest from the 
pulsar (see Figure 6). 

3. DISCUSSION 

In Figure 6, we illustrate the relative contributions of 

the thermal and non-thermal components as a function of 
radius. The thermal shell is clearly present. Note, how- 
ever, that this thermal emission represents less than 4% 
of the observed (absorbed) flux of the remnant in the 
0.5 — 10 keV band, explaining why it has proven so elu- 
sive to date. Figure 7 displays the radial dependence of 
the shell temperature and Neon abundance in the rem- 
nant shell both are flat. In contrast, the power law 
spectral index F monotonically increases from F = 1.9 at 
the remnant center to F = 2.8 at the edge (and F = 2.9 
in the eastern lobe). This latter result is generally consis- 
tent with earlier work (Torn et al 2000; Slane et al. 2002) 
although, owing to the much larger number of photons 



collected and the spectral-spatial decomposition we have 
performed, the measurements extend to a greater radius 
and differ in detail, remaining flat between 2.2' and 3.2'. 

The widths of the annuli provide largely independent 
measurements of the spectrum at each radius given the 
size of the PSF. Some contamination is present from the 
inner, brighter rings to the outer ones, however. Since 
the dominant power law is flatter in the central regions, 
photons scattered to the outer regions bias the power 
law measured there. We have determined the degree of 
contamination from power law photons in each region by 
using the PSF binned in 0.5' annuli and the geometrical 
intersection of those annuli with each extraction annulus 
(the minimal dependence of the PSF on photon energy 
allows us to ignore this complicating factor). The cor- 
rections to the measured power law slopes reported in 
Table 2 are modest, albeit outside the formal errors. In 
the central 0.5' circle, the power law spectral index falls 
from 2.02 to 2.00, while in the first four annuli beyond 
this, the indices changes as follows: 2.29 becomes 2.38, 
2.45 becomes 2.50, 2.55 becomes 2.58, and 2.73 becomes 
2.78. Thus, the rise in spectral slope is slightly steeper 
than shown in Figure 7. PSF contamination effects are 
negligible for the thermal component since the count rate 
is so low. 

The parameters we derive for the shell are in excel- 
lent agreement with those inferred in the early work of 
Bocchino et al. (2001) and with those found from a por- 
tion of the shell in the deep Chandra observation of 3C58 
(Slane et al. 2004): the temperature and Neon abun- 
dance are identical at kT = 0.23 keV and 3.2 solar, re- 
spectively. While these previous authors have cited the 
Neon overabundance as evidence that we are seeing emis- 
sion from the ejecta, a growing body of work (Cunha et 
al. 2006 and references therein) suggests that the solar 
abundance of Neon has been systematically underesti- 
mated by a factor of 2.5 to 3. In that our fits do not 
require the statistically marginal enhancement of a fac- 
tor of two in the Mg abundance reported by Slane et 
al. (2004), we argue that the provenance of the thermal 
emission remains an open question. 

The total observed X-ray luminosity we derive for the 
shell component is 5.9 x 10'^^(i| 2 ^^'g^ ^^^i within the 
range found in Bocchino et al. (2001). This is signif- 
icantly below the upper limit derived by Seward et al. 
(2006) for the missing Crab Nebula shell. The total mass 
is ^ O.77M0 (for a mean molecular weight of 0.6 which 
may not be appropriate if the bulk of the radiating ma- 
terial is ejecta). 

Evidence has been accumulating for some time that, 
despite the apparently robust conclusions drawn from the 
historical records that 3C58 is coincident with SN1181 
(Stephenson and Green 1999), the remnant properties 
are inconsistent with such a young age. Bietenholz (2006) 
provides a comprehensive summary of the arguments 
against an association, as well as a list of alternative sce- 
narios. Most of the evidence suggests a remnant age of 
~ 3000-5000 years (e.g., Chevalier 2005). Our spatially 
resolved image of the shell adds further to the case for 
an age > 820 yrs. 

Using the X-ray emitting mass we have derived and the 
analytic models of PWN evolution from Chevalier (2005) 
we can set limits on the apparent age of the remnant. 
Prom Eq. 27 of Chevalier (2005), we can combine the 
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Fig. 8. — X-ray image of the pulsar wind nebula and the thermal shell in the supernova remnant 3C58. These images are the mosaic of 8 
observations acquired with the XMM-Newton EPIC cameras. Left) the 1.0 — 2.0 keV exposure-corrected image with contours illustrates 
the morphology of the pulsar wind nebula and "lobe" regions. This image is scaled to enhance the fainter nebula emission with contours 
overlayed. Right) The same data but in the 0.5 — 1.0 keV energy band with the non-thermal contribution subtracted (see text). The 
intensity is scaled linearly and the seven contours are equally spaced. The symmetric, filled shell-like structure is ~ 5.'6 in diameter. 



observed supernova remnant radius of 2.8 d^,2 pc and 
the pulsar timing-derived parameter E — 2.7 x 10^^ erg 
s~^ to obtain the age in units of 1000 yrs, 

^3 « 1.8 E-^-^ 

where M^j is the ejected mass in units of solar mass Mq. 
If we make the extreme assumption that all of the ob- 
served Q.ITMq of X-ray-emitting gas is ejecta, and that 
it represents some fraction < 1 of the total ejected mass 
of the supernova, we find ^3 ^ 1.62 E^^''^. Thus, the pul- 
sar must be at least twice the age of SN1181 for reason- 
able explosion energies and, for a more plausible ejected 
mass of ~ 5Mq, the age is 3400 yrs. Alternatively, if 
we assume all of the observed gas has been swept up by 
the expanding PWN, from Chevalier's (2005) eq. 28, we 
find t — 5100 yrs (cf. the pulsar's characteristic age of 
5390 yrs). 

The pulsar's location with respect to the circular shell 
of thermal X-ray emission provides a potential test that 
could resolve the age controversy. The location of the 
pulsar is marked with a crosshair in Figure 8. It is ap- 
parent that this is not the center of the shell; a circle 
best matching the observed thermal emission is centered 
27"±5" west (and slightly north) of the pulsar's location. 



Bietenholz. M. F. 2006, ArXiv Astrophysics e-prints, 

arXiv:astro-ph/0603197 
Bocchmo, t vvarwicK, n. S., Marty, P., Lumb, D., Becker, W., & 

Pigot, C. 2001, A&A, 369, 1078 
ChevaUer, R. A. 2005, ApJ, 619, 839 

Cunha, K.. Hubenv. I.. & Lanz. T. 2006, ArXiv Astrophysics e- 

prints, arXiv:astro-ph/0606738 
Fesen, R. A. iyo3, ApJ, 2ru, Lo3 

Gaensler, B. M., Arons, J., Kaspi, V. M., PivovarofT, M. J., Kawai, 

N., & Tamura, K. 2002, ApJ, 569, 878 
Helfand, D. J., Gotthelf, E. V., & Halpern, J. P. 2001, ApJ, 556, 

380 

Jansen, F., et al. 2001, A&A, 365, LI 
Lumb, D. H., et al. 2004, A&A, 420, 853 



For an age of 821 years, a distance of 3.2 kpc, and spher- 
ical expansion (an assumption supported by the fact we 
see no thermal emission in the eastern elongation of the 
remnant), this requires a two-dimensional pulsar veloc- 
ity of 500(i3.2 km s~^. While not impossible for a young 
neutron star, it is worth noting that if the object has the 
same tranverse velocity as the Crab pulsar, the remnant 
age would be 3750^3.2 yrs, consistent with other esti- 
mates. We further note that the implied trajectory for 
the pulsar aligns to within a few degrees of the long axis 
of the nebula, behavior commonly seen among PWNe. If 
Chandra survives for another decade, it will be possible 
to resolve this issue directly, as the implied proper mo- 
tion of 0.5" for the high- velocity scenario will be directly 
measurable. 



This research is supported by NASA LTSA grant 
NAG 5-8063 to EVG and by grant SAO GO3-4026B to 
DJH. This research has also made use of data obtained 
from the High Energy Astrophysics Science Archive Re- 
search Center (HEASARC), provided by NASA's God- 
dard Space Flight Center. 



Reynolds, S. P. & AUer, H. D. 1988, ApJ, 327, 845 
Saxton, R. D. 2003, XMM-SOC-CAL-TN-0023 
Seward, F. D., Gorenstein, P., & Smith, R. K. 2006, ApJ, 636, 873 
Slane, P. O., Helfand, D. J., & Murray, S. S. 2002, ApJ, 571, L45 
Slane, P., Helfand, D. J., van der Swaluw, E., & Murray, S. S. 2004, 
ApJ, 616, 403 

Stephenson, F. R. , & Green, D. A. 1999, Astronomy and 

Geophysics, 40, 27 
Torii, K., Slane, P. O., Kinugasa, K., Hashimotodani, K., & 

Tsunemi, H. 2000, PASJ, 52, 875 
Turner, M. J. L., Briel, U. G., Ferrando, P., Griffiths, R. G., & 

Villa, G. E. 2003, SPIE, 4851, 169 



Thermal X-ray Shell Around SNR 3C58 



7 



TABLE 2 
Spectral Fits and Fluxes 



Model 








Annulus 








Parameter 


(oro-o.'5) 


(o.'5-iro) 


(I'.O-l'.b) 


{l.'5-2ro) 


(2.'0-2r5) 


{2'.5-3'.0) 


(3r0-3.'5) 



Nn (1021 
F (spectral index) 
PL Flux^ pn 
PL Flux* MOS 
x'(DoF) 



kT (keV) 

Ne ([Ne/Ne©]) 

kT Flux* pn 

kT Flux* MOS 

x'(DoF) 

Total Flux* pn 

Total Flux* MOS 



4.16(4.09-4.24) 
2.02(2.01-2.04) 
3.43 X 10-12 
4.14 X 10-12 
404.01(442) 



0.24(fixed) 
3.0(fixed) 



399.37(443) 
3.43 X 10-12 
4.14 X 10-12 



4.16(fixed) 
2.29(2.27-2.31) 
2.80 X 10-12 
3.15 X 10-12 
207.01(251) 



0.24(0.20-0.29) 
5.5(3.0-9.4) 
0.27 X 10-13 
0.18 X 10-13 
373.58(414) 
2.82 X 10-12 
3.20 X 10 



4.16(fixed) 
2.45(2.43-2.48) 
2.06 X 10-12 
2.38 X 10-12 
138.71(196) 



0.22(0.20-0.24) 
4.3(3.6-4.8) 
0.63 X 10-13 
0.57 X 10-13 
294.33(365) 
2.12 X 10-12 
2.41 X 10-12 



4.16(fixed) 
2.55(2.51-2.59) 
1.50 X 10-12 
1.67 X 10-12 
129.18(152) 



0.22(0.21-0.23) 
3.1(2.7-3.4) 
1.16 X 10-13 
1.16 X 10-13 
345.91(319) 
1.62 X 10-12 
1.78 X 10-12 



4.16(fixed) 
2.73(2.67-2.78) 
8.40 X 10-13 
9.60 X 10-13 
80.28 (101) 



0.23(0.22-0.23) 
3.0(2.8-3.3) 
1.39 X 10-13 
1.63 X 10-13 
324.62(252) 
9.76 X 10-13 
1.12 X 10-12 



4.16(fixcd) 
2.81(2.70-2.92) 
3.86 X 10-13 
4.40 X 10-13 
36.98(66) 



0.23(0.22-0.24) 
3.6(3.2-4.0) 
0.96 X 10" 
1.05 X 10-13 
140.76(164) 
4.79 X 10-13 
5.40 X 10-13 



4.16(fixed) 
2.93(2.75-3.12) 
2.05 X 10-13 
2.48 X 10-13 
39.66(52) 



0.24(0.21-0.28) 
3.1(2.2-4.4) 
0.26 X 10-13 
0.26 X 10-13 
117.61(125) 
2.34 X 10-13 
2.62 X 10-13 
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Note. — Uncertainties are 90% confidence for two interesting parameters. 
* Absorbed flux in tlie 0.5-10 keV band in units of (ergs cm"^ s— "'^). 



